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Kelly type joint measurement model for qubits. We investigate the effect of initial state of the 
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the system. Two approaches for characterising the quality of unsharp measurement and the resulting 
measurement uncertainty relations are considered. The corresponding measures of unsharpness are 
connected for the case where both the measurements are equally unsharp. The connection between 
the POVM elements and symmetries of the underlying Hamiltonian of the measurement interaction 
is made explicit and used to perform joint measurement in arbitrary directions. Finally in the case 
of three observables we derive a necessary conditon for the approximate joint measurement and use 
it show the relative freedom available when the observables are non-orthogonal. 
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I. INTRODUCTION 



Quantum mechanics does not allow the joint measurement of non-commuting observables. The uncertainty principle, 
usually described as lower bound on the product of standard deviations of the outcome statistics does not really capture 
this complementary feature as it (uncertainty principle) is a statement about the quantum state of the system and 
does not relate to an actual situation involving apparatus which can attempt a joint measurement. 

As intuition suggests one has to allow for some degree of imprecision in order to make room for a notion of joint 
measurement of non-commuting observables. In quantum mechanics of a single system observables correspond to 
self-adjoint operators with the outcome statistics given by the Born rule , i.e if Pi is the projector onto the eigenspace 
corresponding to the eigenvalue then in a measurement of the observable the probability of obtaining the ith outcome 
when the system is in state p is Tr(pPi). However this is not adequate. For example, suppose a system is allowed to 
interact with another system (ancilla) for a while and then measurement is done on the combined system. To describe 
the probabilities of obtaining various outcomes for such a measurement for different states of the original system 
before the interaction one has to replace the projectors above by positive operators £7, (actiing on the states of the 
original sysytem) with £\ = I for normalisation. Such measurements are called POVM (positive operator valued 
measures) measurements. The most general measurements possible in quantum mechanics can be described through 
POVMs. It turns out that the idea of imprecise or unsharp measurements can be appropriately investigated if one 
considers instead of projective measurements, POVMs. The observables corresponding to projective measurements 
form a subclass of those described by POVMs and are called sharp. 

The first model for approximate joint measurement of position and momentum was given by Arthur and Kelly in 
(1965) by generalising von Neumann's model for measurement. ([2]) In the Von Neumann's model for measurement, 
the position observable of the object is measured by coupling it to the momentum P p of a probe system via the 
interaction evolution U = e~i x Q® Pp and using the position Q p of the probe as the readout observable. The idea of 
Arthur-Kelly was to couple two such probe systems respectively to the position and momentum of the system and 
then perform measurement on the commuting meter observables of the probe to gain information about the position 
and momenta of the system. This is in fact an unsharp joint measurement and inspired the later development of the 
formalism . EPR-Bcll argument for unsharp realities was given in [17J ; the bounds on joint measurement of qubit 
observables were derived starting from the condition of operational locality in [S] ; the formalism for qubit observables 
and relevant inequalities for the two observable case were developed in [7], [H],[E], [ID]- The connection of approximate 
joint measurement with Bell inequalities were investigated in |1 9 j : joint measurements by quantum cloning has been 
invesigated for example in [20] and [5T].It has also been used to estimate the expectation values of qubit observables 
in |16j ; joint unsharp measurement of position and momentum was treated in [?j and [5]. 
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In this paper we consider an Arthur-Kelly like model for approximate joint measurement separately, of two and 
three non-commuting qubit observables . The probe systems we employ are continuous as typically in experiments like 
Stern-Gerlach deviations in position or momenta of particles carrying spin are measured to obtain information about 
the system spin state. We consider characterisation of unsharpnesss or quality of approximation by two ways existing 
in literature. Firstly by considering closeness of the marginal probability distribution of the probability distribution 
of joint measurement to that of the sharp observable being approximately measured .Secondly by considering suitably 
defined closeness of the observables themselves with the meter observables in the Heisenberg picture. We show that for 
a symmetric joint measurement where the two marginal probability distributions are equally close to the corresponding 
sharp probability distributions, the two measures of closeness are proportional. Error-disturbance relationship does not 
seem to hold for the measure based on the Heisenberg picture for our choice of the pointer observable. We also prove a 
lemma showing the connection between joint measurement ,and the symmetries of the underlying Hamiltonian of the 
measurement interaction together with that of the initial detector states. This is then used to perform approximate 
joint measurement in arbitrary directions. Moving on to the case of three-observable joint unsharp measurement we 
prove a simple necessary condition that is sufficient for the case of three orthogonal observables. The case of three 
orthogonal observables also generalises the validity of the known necessary-sufficient condition for this case. This 
condition is derived from certain geometrical considerations based on the so called Fermat-Toricelli point. We show 
that it yields the known two-observable bounds in the limit that the measurement of one of the observables is pure 
guessing of the value of the observable. Numerical analysis is performed for the two observable case to show the 
validity of measurement uncertainty relations, transition from POVM to projection-valued measurement and also the 
effect of the joint measurement on the system. Similar analysis is also performed for the three-observable case. 

The paper is organised as follows . In section II we introduce the joint unsharp spin measurement and the original 
Arthur-Kelly model given for the joint measurement of position and momentum. Measures characterising the quality 
of approximations to the sharp observables being approximated are introduced in section III . The corresponding 
measurement uncertainty relations are discussed. In section IV we introduce the Arthur-Kelly like model for qubits 
that we have considered and derive the final state of the sysytem and meters after the measurement interaction. 
Approximate joint measurement of a x and a y is considered in section V. Section VI deals with the effect of the initial 
detector states on the joint measurement and post-measurement state of the system. In section VII we explore through 
a lemma the connection between the symmetries of the underlying Hamiltonian of the measurement interaction and 
initial detector states on the joint measurement. This is used in section VIII to do joint measurement of spin in 
arbitrary directions. The corresponding POVM elements are also calculated and matched with the orthogonal case in 
section V. In section IX we compute the spin direction fidelities, which were introduced as a measure of quality of 
approximation in section III, for the model considered. It is compared with another measure based on distance between 
outcome probabilities of sharp observables and their unsharp approximations. Approximate joint measurement of 
three qubit observables is considered in section X. A necessary condition on the parameters of the marginal POVM 
elements is derived by geometric considerations involving the Fermat-Toricelli point. The suffciency of the condition 
is explored and the restrictions placed by it investigated in the context of joint measurement through an extension of 
the model considered in section IV to three detectors. 



II. AN INTRODUCTION TO JOINT UNSHARP SPIN MEASUREMENTS AND ARTHUR-KELLY 

MODEL 

Suppose we want to measure the spin of a spin-| particle along a direction given by the unit vector h .If the 
system is described by the density matrix p, then probabilities for obtaining outcomes ± are given in the POVM 
formalism, respectively by Tr(pA + ) and Tr(pA_). The POVM elements A + and A_ are called effects. One says 
that the observable A is charaterised by specification of the map Wi —> A4 where Ui belongs to the set of outcomes 
with i = +,— (e.g here u + = + , w_ = — ) . If A + = |n, +){n, +\ i.e a projector into the up state of a.h then the 
measurement is called a sharp measurement of a.ri and the observable characterised by A + (A- = I — A + ) a sharp 
observable. Else the measurement and the corresponding observable is called unsharp. One way such a measurement 
can occur is the following. Suppose some one wants to do a standard (sharp)measurement of a.h but because of some 
error in his setup there is a finite probabilty of registering the outcome - when the system is actually in the state 
|+) and vice versa. This will be an unsharp measurement with the probabilities of error given by Tr(A— and 
Tr(A + \— )(— |). Such a situation can arise if for example in a Stern-Gerlach setup , the beam passing through the 
inhomogeneous magnetic field is poorly collimated initially. [3J. 
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A. Joint measurability 

Two observables are said to be jointly measurable if there is a measurement scheme that allows the determination of 
values of both the observables. This means that the POVM describing the measurement scheme contains the POVM 
elements of the two observables as marginals. In this way it is ensured that there is a joint probability distribution 
corresponding to the pairs of different values of the two observables for each state. Two observables Ti and Y2 are 
jointly measurable if there is an observable G : u>ij —¥ G^, i,j — ±, such that, 

T + = G ++ + G + - 
= G ++ + G_ + 
ti = G_+ + G— 

Ti = G+_ +G__ (1) 

with the POVM elements f £(for j = 1, 2) satisfying , f \ + T 1 = / and f \ + f 2 _ = I . The outcomes uj v] of G 
can be taken to be the pairs (aij, Wj). The map implies that the probability of registering outcome , say (+,+) is given 
by Tr(pG++) and so on. The extension to three observables is done in the obvious way. 

Two observables Ti and T2 are said to commute if T; Y~ = Y^Y* Vi, j = ±. The possibility of joint measurability 
is guaranteed by the following theorem which we state without proof. ([H]) 

Theorem 1: A pair of sharp observables is jointly measurable iff they commute . Commutativity of unsharp 
observables is sufficient but not necessary for joint measurability. 



B. Arthur-Kelly Model 



As mentioned earlier Arthur and Kelly gave a model for joint measurement of position and momentum. 2J In this 
model a quantum object is coupled with two probe systems which are then individually measured to obtain information 
about objects position and momentum. The coupling to probes is based on von-Neumann's model of measurement. 
They showed that this constitutes a simultaneous measurement of position and momentum in the sense that the 
output statistics reproduce the expectation values of the object's position and momentum. 

The position Q and momentum P of the object are coupled with the position Qi and the momentum P 2 of the 
two probe systems respectively which serve as the readout observables. Neglecting the free evolutions of the three 
systems (assuming that the measurement interaction dominates during the short time in which it acts) the combined 
time evolution is described by , 

The coupling constants A and k represent the interaction strength and can be absorbed into a rescaling of the pointer 
observables Qi and P%. If is an arbitrary input state of the object and l^i) and |\l/ 2 ) are fixed initial state of the 
probes (well-behaved and zero expectations for each of the probe's position and momentum ) , the probabilities for 
values of Qi and P2 to lie in intervals of length Xx and ny are determined through a joint observable G T [4] , analogous 
to the one defined in eqn. ([l]) , through 

(4>\G T (XxY)\^) := (i>\ ® ® (V 2 \UHI® Qi (Ax) ® P 2 {Ky))U\f) ® |* 2 ). (3) 

The marginals are given by Gj(X) := J2 y ey G T (XxY) and G^iY) — J2xex G T (XxY) , where X and y represent 

the set of all possible values respectively of Qi and P 2 -They play the same role as T 1 and T 2 in[l] , with for example 
the probability to obtain Q\ in an interval Ax (about x = 0)after the approximate joint measurement given by 
(4>\Gf (X)\ip\) with X = [— 4p, 4p] . The cost of the joint measurement is an increase in the variance of the marginals 
given by, 

A(Gf,^)A(G 2 n ,V)>^ (4) 



where, A(G?» = \[w%SI>) ~ 
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III. QUALITY OF JOINT MEASUREMENT AND MEASUREMENT UNCERTAINTY RELATIONS 

The notion of approximate joint measurement naturally demands a measure of proximity to the sharp observables 
being approximately measured. The restriction on the measures corresponding to non-commuting observables leads 
one to measurement uncertainty relations. In this paper we use two such approaches to characterise proximity. 

A. Closeness of probabilities 

Proximity between two observables can be characterised by the distance between the corresponding probability 
distributions for all states. Thus distance between two observables A and B can be defined as 

D(A, B) := max,-sup r |ir[TA,-] - trlTBj]] (5) 

where Aj(orBj) corresponds to the POVM element (or effect) of the observable A (or B) associated with the measure- 
ment outcome j. T is the density matrix of the sysytem.JTj 
For two single-qubit observables T°^ a } T l3:b with, 

t«> s =±(aI + aJ) (6) 

T^= l -{pi + b.ff) (7) 

with (a, a) and (/3, b) G 1Z . Now for observables A , B with respective set of effects {A + , A-} , {£>+, B^} we have 
\tr[TA+] - tr[TB+]\ = \tr[TAJ\ - tr[TB-]\. Taking T to be the '+' state or '-' state of a. .^Zrln according as 
(a — 0) > o or (a — fi) < o one has, 

D{T«>^r^) = \\\a-b\\ + \\a-p\ (8) 

This shows that the distance of a certain unsharp observable T Q > a from any sharp observable T 1,n is minimum when 
a is along n. 

1. Unbiased Observables 

Observables of the form f x - 3 are called unbiased. As f ^ s = ^M^I + w) + ^^-^(I - therefore 
the probability of occurence of outcomes +(or -) for the initial state | (the maximally mixed state) is given by, 
Tr[^.\{I + ^p] = \ ^or Tr[|.|(J — ^pfp)] = 5) • ^° t nese t w0 probabilities are same. Hence the name 'unbiased'. 

For such an observable, both the outcomes are equally likely for a maximally mixed state . Also, the expectation 
value of the unsharp measurement when the system is in state T is given by (a. a ) u = l.tr{T\{I + a.a)) -l.tr (T\(I- 

a.a)) = \\a\\(<j.a) . Again as, D(T 1,a , T 1 ' ) = |(1 — ||a||), ||a|| itself serves as a measure of proximity . We will often 
approximate sharp observables with unbiased unsharp ones. 

2. Measurement uncertainties 

We will choose jointly measurable observable pairs (T a '°, T' 9 ' 6 ) to approximate the sharp pair (T 1 '^, T 1 '*). The 
necessary and sufficient conditions on a, ft, a, b so that the first pair is jointly measurable in the sense of[l]is formulated 
in [S],[S], [ID]- When the observables are unbiased the conditions simplify to, 

||3+6j| + ||3-6||<2 (9) 

It was further shown in [7] that, 

D(T a ' 3 , f x,ft ) + D(t p,i , f x '™) > 2D Q (10) 

with, Do = -^(cos(6/2) + sin(9/2) — 1) , 6 being the angle between a and b . The conditions for attainment 2Dq 
were also spelt out . The approximate observables should be unbiased, i.e a = P = 1 in this later case. The other 
conditions imply that for optimality, a and b lie along h and m respectively, only when the latter are orthogonal [7]. 
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B. Closeness of observables 



A completely different approach was taken in references [T| and [IT] in the context of the original Arthur-Kelly 
model to give a formulation of the complementary nature of the approximate joint measurement process. If p x and jj, p 
denote the pointer observables used to measure the system position and momentum respectively , then the retrodictive 
error operators are defined as , 



the predictive error operators as, 
and the disturbance operators as, 



ex, = tfxf — Xi, £Pi = fip f - Pi (11) 
£Xf=fixf- Xf , ep f = fi p f - p f (12) 
$x = x f - Xi, e P = p f - pi (13) 



where operators Of, appearing on right stands for the final Heisengberg picture operator O after the measurement 
interaction U, i.e, Of — U'OiU where Oj is the Heisenberg picture operator at the moment the interaction starts. 
Various errors were then defined by taking the square root of expectation of the square of the operators defined above 
and taking supremum over the system states. For example the maximal error of retrodiction was defined as , 

A e « = sup WeH ^ t ({$®V 1 ®* 3 \eZ i 2 \i/>®V 1 ®9 a ))*, (14) 

and similarly for A e i p , A,,^, A e ^ , Ad x and Ad p (see [1] , [11] for details) .Measurement uncertainty or 'error 
principle' was shown to hold in the form of, 

h 

A e i x A e i p > — , (15) 

h h 

p — 2' ^eix^dp > ^' 

&ef x A dp > | (17) 

and extensions of the above uncertainties in the obvious way. One of the important features of this approach is the 
difference between error of retrodiction and that of prediction .It was shown that these are not the same as long as 
there is a finite disturbance. 



A efx A eft) > -, A ax A dp > -, (16) 



C. Qubit Observables 



In a similar spirit for the case of approximate joint measurement of qubit observables through an Arthur-Kelly like 
process fidelities were defined in the Heisenberg picture that would provide a notion of direction of spin |12) of the 
system. As in the previous case distinction was made between errors of retrodiction and prediction. In this paper we 
consider only the type 1 measurements considered by the author. We next consider the fidelities as defined by the 
authors. The retrodictive fidelity is defined as, 

Vi - tnf lx)eHsys {(^\ <g> + Si.ri f )\il>) |x» (18) 

where, Si = S and Hi = h are the initial values of the Heisenberg spin and pointer observables and Sf — U*{I ® 
Si)U and rif = (rii (8> 1 S )U be the final Heisenberg spin and pointer direction observables after the measurement 
interaction. 

The predictive fidelity is defined as, 

Vf = in f\x)eH By3 M ® (x\\(nf-S f + Sf.ti f )\i>) ® |x»- (19) 
The measurement disturbance by, 



m = inf lx)eHsys m ® (xl^SfA + Si.SfM ® |x)). 



(20) 
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The intuition behind the above definitions is classical in the sense that it considers alignment of initial or final 
spin vector and initial or final pointer direction. But the above fidelities were used to define maximal rms error of 
rctrodiction, 



a 2 2 

A ei S = sup lx)eHsys ((ip\ <g) (xWVirif ~ Si\ ® |x)) 

= (.s + s 2 -r^, (21) 

maximal rms error of prediction, 

- 2 5 

A ef S = sup lx)eHsys ((^\ ® (xWVfftf ~ S f \ |V>) ® |x}) 

= + (22) 

and maximal rms disturbance 

A d S = S up lx)eHs J(iP\ ® (x||S/ - S/|V) ® | X )) 5 

= N /2( s + s 2 -r^, (23) 

where the spin s = \ for our case and \r)irif — Si\ 2 = {rjiUf — Si) 2 . The quantities A e iS , A e fS, AdS were expected 
to play the same role for these measurements as similar quantities defined for the original Arthur-Kelly model. The 
authors in ref.[T2] also showed that the retodictive and predictive fidelities reach their optimal values of s — 0.5 when 
the Krauss operators of the measurement POVM have a spin-coherent form defined using spin coherent states. A 
physical model for such a measurement was given in |13j . 

No measurement uncertainties were derived in |12j and the question was left open for further investigation. 

IV. THE MODEL 

We consider an instantaneous coupling interaction with the help of the Hamiltonian of the form, 

H = -(q 1 <Z>(T x + q 2 (S)<j y )5t (24) 

(in h=l units) Possible coupling constants in the above equation like A and k in equation [2] have been absorbed 
by rescaling q\ and q 2 . Like the original Arthur-Kelly interaction ^ the idea is to entangle the detectors with the 
system through H and then perform a projective measurement of pi , p 2 to obtain the spin information. Now, as 
a consequence of the Ehrcnfest theorem the average momentum change of a particle carrying spin that experiences 
the above interaction is given by < p\ >~< a x > and < p 2 >=< a y > . Thus, an ensemble of particles whose spin 
state is | + x > and which has a symmetric distribution of p\ before the interaction will have a greater probability 



of having a positive p\ after it. The signs in equation (24 1 have been chosen so that this is true for both x and y 



directions. The signs thus allow us to map the four quadrants of the momentum plane (pi,p 2 ) to the four outcomes 
of joint measurement and take the signs of the momenta to correspond to the outcomes (+,+) , (+,-),(-,+) ,(-,-) of 



the joint measurement. Note that for a Stern-Gerlach situation the two terms in eqn. (24) should have opposite signs 
to satisfy divergenceless of magnetic field. 

Models similar to above have been considered before for example in [3] , [2] , [TS] ■ As shown in [H] this model 
naturally arises in the context of a Stern-Gerlach experiment with a linear magnetic field. 



We further assume that the measurement interaction (24 1 is strong enough to dominate the other parts of the 



Hamiltonian during its presence (e.g, the kinetic energy part ). In the Stern-Gerlach context this would mean to 
assume the atoms carrying spin to be sufficiently massive. 



The unitary evolution corresponding to the Hamiltonian of equation ( 24 1 is given by, 



U = exp(i(qi ® a x + q 2 ® a v )) (25) 
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On direct expansion, U turns out to have the simple form, 

U = e(gi, g 2 ) ® 1* ®o-a:+ fl(9l,92) ®0"d 
with, e(qi,q 2 ) = cos{^J (qf + qfj) 

/(9i,92) = *9i 7=3 ^ 

x/Cgx + 9|) 

5(91,92) = «92 * „ ? . (26) 

V(9?+9 2 2 ) 

Thus, the final state of the system is given by 



\^f) = L ,,=-00 e (9l> 92)191,92)^1 {qi)^2{q2)dqidq 2 ® \ X ) 



J 9l,92=-oo 
p+00 

f +00 



+ f qi ™ 2 =-oc 92)|9l, 92)^1 (9l) "02 (92)^91^92 ® C^lx) 



+ X 91 ^ 2 =-oo 5(9l,92)|9l, 92)^1(91)^2(92)^91^92 ®CT„|x), ( 27 ) 

with initial state being, jt/;,) = ® |i/> 2 ) ® Let, p = l^/XV'/l • 

V. APPROXIMATE JOINT MEASUREMENT IN ORTHOGONAL DIRECTIONS 

In this section we consider the joint measurement of a x and a y .We choose the observables p\ and p 2 to serve as 
meters. As mentioned before, (pi > 0,p 2 > 0) is taken to correspond to the outcome (a x — l,<r y = 1) = (+, +) of 
joint measurement, (pi >0,p 2 < 0) to (+,-), (pi <0,p 2 > 0) to (-,+) and (pi <0,p 2 < 0) to (-,-). 

After the interaction U between the system and meters , projective measurement is performed separately on the 
meters p\ and p 2 . The probability of obtaining outcome (pi,p 2 ) is given by, 

p(Pi,P2) = Tr{\p 1 ,p 2 ){p u p 2 \Tr s {p)) 

= (|e°| 2 + |/T + |5°| 2 )(X|X> + 2Re(f°e°*)( X \a x \ X )+ 

2Re(g°e°*)( X \a y \ X ) - 2Im(g° f*)( X \a z \ X ) (28) 

with e°, / and g° representing respectively the fourier transforms of etpiip 2 , f^i^i &ud gipiip 2 - 
The initial states are taken to be Gaussian, 

fctaO = (29) 

= (^=e-^)i (30) 

We have chosen the initial states to be even in qi , q 2 . Now as fourier transform of a purely real even function is a 
purely real even function and that of a purely imaginary odd function is a purely real odd function we have / odd 
in pi and even in p 2 , g°the other way around and e° even in both Also, each of them is real. Thus, the a z term in 



(281 vanishes. Henceforth in this paper we refer to these properties of e , / and g as "parity properties" 



We have for the probability of outcome (+,+) from equation (28), 

P(P! > 0, P2 > 0) = J™ =0 / p 7 =0 {(|e°| 2 + |/T + |5°| 2 )(x|x) + 2ite(/°e°*)< X |<7 x | X > 

+2Re(g e *)( X \a y \ X )}d Pl dp 2 . (31) 

One also of course has to satisfy, 

p(px>0)+p(pi<Q) = l, 

which yields due to the "parity properties" , 

f°° 1 

/ (\e°\ 2 + \f \ 2 + \g \ 2 )d Pl dp 2 = -. (32) 



From equation (32) and "parity properties" we have, 



p( P i > 0,p 2 >0) = ( x \ I - + ^ + b ^L\ x ) = ( X \G ++ \x) 



with, Gj 



(| + ^ + and 



+00 p-\-oc 



4(/ u e> Pl dp 2 , 



Px=Q Jp 2 — — oo 
oc p4-oo 



b' = 



/ / 4(g°e )d Pl d P2 . 



(33) 

(34) 
(35) 



For the other outcomes, we have from consideration of the corresponding momentum probabilities 

G, = 



G- 



I 


a 


'cr x 


b'a y 


4 


+ - 


4 


4 


/ 


a 


'o x 


b^y 


4 




4 


4 


/ 


a 


Ox 




4 




4 


4 



(36) 



From equations (33) and (36) we get the marginal unsharp observables as, ( see eqn. (|I 



T| = \{I±b'a y ) 



(37) 



Thus, we see that the approximate observables T 1 , T 2 are unbiased and a' and b' themselves serve as measures of 
proximity to the sharp observables ^(I ±a x ) and ^(I±a y ) , respectively. 



In section VII we will see how the approximate joint measurement, characterised by the marginals in eqn. (37) 



arise as a consequence of symmetries of the Hamiltonian in eqn, (24) and that of the initial detector states "01 and 1P2 



( given in eqns. (29) and (30) 



A. Post-measurement state of the system 

In this section we consider the effect of approximate joint measurement of a x and a y on the state of the system. We 
consider the case when the result of the measurement on the meters turns out to be : pi,P2 G [0, 00], i.e the outcome 
(+,+) of the joint measurement has clicked. 

The reduced unnormalised state of the entire system plus the two apparatuses when we concentrate on the pointer 
values Pi,P2 G [0, 00], is given by, 



+00 p+00 



p' 1= Jp' 2 =0 



\p'l,P2)(Pl,P2Wldp' 2 <g)ls \lpf) 



(38) 



Suppose the initial system state is p\ = \{I -\- xa x + ya y + za z ). To obtain the post-measurement statistics for the 
system the detector parts are traced out and we have, after normalisation, the following density matrix of the system: 



Ps 

where, x' 

y' 



-(I + x'a x + y'a y + z'a z ) 

^+4yc fg .+2x(\-2c gg .) 

xa' + yb' + 1 
%+4xc fg .+2y(\-2c fr ) 

xa' + yb' + 1 
2z(|-2( C// . +c flg .)) 
xa' + ybl + 1 



(39) 
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with, 



c fg* 
c ff* 



+ 00 /*+oo 



f°g°*dp[dp' 2 , 



P f i=0 J p' 2 =o 

+ 00 p-\-OQ 



p' 1= Jp' 2 =0 

OG p-\-OQ 



\f\ 2 dpWi- 



Pi=0 Jp' 2 =0 



\g°\ 2 dp[dp' 2 . 



(40) 
(41) 
(42) 



VI. EFFECT OF INITIAL DETECTOR STATES 



The form of the POVM elements in eqns. ( 33 ) , ( 36 1 only depends on the "parity properties " and not on the Gaussian 



form of the initial detector states given by eqns. (29) , (30). For computing the post measurement statistics, the 
Bloch vector of the initial state p\ of the system is chosen to be ( ^ , ^ , 0) ( so that the state uncertainty product 



Ao~ x .Ao~ y is minimum) 



The integrals for computing the coefficients cj g *, ctf and c gff » have been done using the Monte Carlo integrator 
included in GNU Scientific Library. The MISER algorithm has been used which uses stratified random sampling 
for doing importance sampling 23j.The Parseval's theorem can be used to convert the integrals in c//*, c gff » to 
configuration space, e.g, 



c ff , = - I l/l ipi dqidq 2 

' 31,92=0 



A. Results 

As we are considering unbiased approximate joint measurement in orthogonal directions , eqn. (|9| applied on the 
parameters a' and b' of the marginals in eqn. ( 37 ) yield , 



b' 2 < 1. 



(43) 



1. Symmetric case 



In this case the standard deviations a and b of the initial momentum wavefunctions in eqns. (29) and (30) are taken 
to be equal. As there is nothing else that differentiates between the detectors one has a' — b' . This is reflected in 
eqn. (34). Both the observables a x and o~ y are approximated equally well. From eqn. (43) we have, a' < -3= = 0.707 . 

Quality: In figure [l] we see that when a and b are close to zero, which corresponds to the fact that the initial 
momenta p\ , p 2 have large spread we have almost no information about the spin state by looking at the momentum. In 
this case a' and b' are both close to zero and we have trivial marginal effects (|) which represent guessing the value of 
the observable with equal probability for + and -. With increase in a, a' increases and touches .628 at a = b = 0.7. (see 
figure 1(a) ). The graph of a' vrs. a is identical to the one obtained in [14] for the same case. 

Noise: As mentioned earlier, the initial system state is chosen so that Ao~ x Aa y — | . In the post measurement 
state this product almost reaches 1, showing that there is considerable increase in noise in the system(see fig. l(b)[ ). 
The state also becomes mixed with the bloch vector magnitude around 0.5 (see fig. 1(c) I . 



b > a, POVM to PVM transition 



a) Let us choose a — 0.1 . The lhs of eqn. (431 starts off at a low value for b = a and gradually increases and closes 
on the bound (1.0) as b becomes much greater than a. (see figure [2(a) ) b much greater than a reflects the situation 
where the initial momentum wavefunction of apparatus 2 is much sharper than that of appararus 1 . As shown clearly 



by figure 2(a) this marks a transition from a POVM measurement to a projection- valued measurement (PVM) in the 



sense that the unsharp measurement of o~ y becomes almost sharp. The requirement of complementarity is satisfied by 
the fact that the unsharp a x measurement becomes almost trivial. The fact that a' and b 1 depend both on a and b is 
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FIG. 1: Effects of initial pointer states on the joint measurement and post-measurement system state for a = b . a) A plot of 
a' vs. the initial detector state standard deviation a , b) A mean plot of the post measurement state uncertainty product of 
the system with a, c) A mean plot of the magnitude of the bloch vector of the post measurement state of the system with a . 




0.5 1 1,5 2 2.5 3 0.5 1 1.5 2 2,5 3 

bwitha-0.1 b 



FIG. 2: Effects of initial pointer states on the joint measurement and post-measurement system state for b > a = 0.1. a)A plot 
of proximities a' and b' vrs. b for a=0.1. The bottom curve shows a' and the top curve shows &', b)A mean plot of the angle 
made by the bloch vector of the post-measurement state of the system with x-axis, given the outcome (+,+) vs. b 



a reflection of the entanglement between pi and P2 brought about by eqn. (25 1. The POVM to PVM transition does 
not have much effect that differ from fig. [l]on the state uncertainty product and bloch vector magnitude. However 
because in order to examine the post-measurement statistics the (+,+) outcome was considered, the bloch vector 
points almost along y direction (see fig. |2(b)[ ). This has been seen to hold indcpendant of the choice of the initial 
state of the system. 

b) Next we choose a = 1.0. The conclusions are almost similar with the exception that the joint measurement 
uncertainty (i.e the l.h.s of eqn. ^ starts off at a much higher value (see fig. 3). 
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FIG. 3: A plot of a' and b' vrs. b, a=1.0 . Bottom curve shows a' and top curve shows b' 



VII. EFFECT OF THE SYMMETRIES OF THE UNDERLYING HAMILTONIAN ON THE POVM 

ELEMENTS 

It was shown in [3] in the context of a Stern-Gerlach type Hamiltonian how the symmetries of the underlying 
Hamiltonian like the one given in eqn. ( 24 ) can be used to do some particular unsharp measurement. The various 
cases considered there can be seen to follow from the following lemma. 

Consider an Arthur-Kelly like measurement process and after the final measurement on the detectors is made. 

Lemma 1: Let H be an Arthur-Kelly like Hamiltonian of the form H = /(<7i, 92) <£> cr x + g{q\, q 2 ) ® °~y , which has 
a symmetry given by [A B, H] — . Here A and B are unitaries acting respectively on the joint detector space 
and the spin space. Let P\ x \ (pi,f>2) represent the probability of obtaining the momenta values p\ , p 2 if the system is 
initially in the state \x) ■ Further let the initial joint detector state \ip) have the symmetry A so that A\if>) = e^lip). 
Then for the initial system state \\) and a new basis in the joint detector space IPxiPa) = ■^■IPuP^) j we have 



P\ x ){j?i,V2) = Pb\ x )(p[,P 2 ) 



Proof: 



P\x)(Pl,P2) = Tr i,%\P[\pi,Pi)]Tr s {t 



®|x)]e 



with P[|?7)] denoting projector on |»y). 

Now replacing e~ lH by (A t ® B^)e~ iH (A ® B) in 145} and using A\ip) = e*^|t/>) we have 



Let , 



P\ x) {px,P2) = Tr 1 , 2 [P[| P i,P2)]Tr s ((At ® B*)e~ iH {ty) ® B\ X )(x\B*)e m {A® B)] 



(IVXVI ® B|x)0cl-Bt)e iH - e~ lH {PU) ® B| X )])e iH = ^C, ® D j5 



(44) 
(45) 

(46) 
(47) 



where the index j runs over a countable set. 

Cj s and DjS are operators acting respectively on the joint Hilbert space of the two detectors and the system Hilbert 
space. 



Using B^ B = / we have from eqns. (45) , (46) and (47): 

Tr s (e- iH P[\iP) ® \ X ))e iH ) = ^ Tr s (Dj)(A i Cj A) . 



(48) 



So,P lx} (p u p 2 ) = J2 Tr ^[\P^P2}(pi,P2\Tr s (D j )(A^C j A)} 
3 

= Y / (P^P2\^C 1 A\ Pl ,p 2 )Tr s (D J ) 

3 

= J2(P'nP2\ C M,P2)Tr s (D J ) 



(49) 
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Again, 



Pb\ x M>P2) = Tr h2 [P[\p[,p> 2 )]Tr s (e- lH PU)®B\x)y H )] 

= Tri, 3 [P[bi,Pa>]Tr a (£ C j ® £>,-)], (from eqn. @ 

j 

= ^(pi,^|Qbi,^)Tr s p 3 ) 



(50) 
□ 



VIII. APPROXIMATE JOINT MEASUREMENT IN ARBITRARY DIRECTIONS 



Let li (i=l,2) denote the operation of reflection in the detector space about g.;-axis. Consider a Hamiltonian of the 
enire system (i.e two detectors plus the qubit) which satisfies \1\ ® a x ,H] =0 and an initial joint state ip of the two 
detectors jointly, that satisfies Ii\i(>) = \ip) ■ Then eqn. (44) yields, 



So, (x\E(p u p 2 )\x) 
i.e E(p 1 ,p 2 ) 
and integrating over p 2 , [E(j>i), a x ] 

So we can write, E(pi) 

where a{pi) and /3(pi) are real nos. 

Hence we have, E(pi > 0) 



P a x | X )(Pl,-P2)- 
(xWxM(jPl,-p2)<T x \x) 

a x E(pi, -p 2 )a x 
0. 

l(a(( Pl )I + (3(( Pl )a x ) 



-( a 'l + p'a x ) 



(51) 



with constant a' and /?'. 

As the Hamiltonian H of eqn. (24| satisfies \I\ ® a x , H ] = and as the initial joint detector state \tp) — \tpi) ® \ip 2 ), 
where and \ip 2 ) are given respectively by eqns. ( 29 1 and ( 30 ) satisfies both Ii\ip) = \ip) and /2IV') = \4>) , it is 
clear from eqn. 51 that how the appproximate joint measurement of a x and a y with the marginals given by eqn. (37) 



arises out of the measurement of p\ and p 2 . 

Further consider an H which in addition to having the symmetry mentioned in the beginning of the present section, 
has a further rotatational symmetry: [R(8) <g> S z (9), H] — 0, where R{9) and S z {9) respectively denote rotation by an 
angle in the detector Hilbert space (i.e, the Hilbert space operation corresponding to rotation in the qi, q 2 plane) 
and the spin space (about z-axis). From eqn. (44) we have, 



P\x)(Pi>P2) = p s z (8)\ x )(p'i,p' 2 ) 



(52) 



with oi,p' 2 ) T = R(d)( Pl , P2 y 



Integrating both sides over the region p\ > 0, p 2 € (—00, +00) we have of eqn. (52), 



P\x)(Pi >0) = PsAxM >0), 
i.e,P s t |x) (pi>0) = P| x> (pi>0). 



Now, using eqn. 51 we have from eqn. 53 



P\xM >0) 
i.e, (xW, > 0)x) 



(x\S z -(a'I + P'<j x )Si\x), 
{X\\{a'l + p'a.h)\x) 



(53) 

(54) 
(55) 



where h = xcos(9) + ysin(9), with 9 being the polar angle of a point in the x-y plane. In the last but one line we have 
also used S z (9) — e _l ~s~ . 



The Hamiltonian in eqn. (24) satisfies both the above mentioned reflection as well as the rotational symmetry 
properties . Choosing the detector state to be a symmetric Gaussian ,i.e, b — a in eqn. (30) wc have the required 
rotational invariance for all angles. Hence, we get an approximate joint measurement of spin in any direction by 
measuring the detector momentum in that direction. 
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A. POVM elements 

Consider approximate joint measurement of <j x and a.n. We have already shown that the marginal probabilities of 
the joint measurement will be given by, 

Pip'x >0) - (x\\{l + a'o x )\x) (56) 
P (p' 2 >0) = ~(l + a'ff.n) (57) 

with p[ and p 2 denoting momenta in x and n direction respectively. (x is along positive x-axis and n along xcos(8) + 
ysin(8) in the momentum plane (pi,P2)) 

1. Angular dependance of e° , f° and g° 

The rotational invariance of the initial state allows one to extract the polar angular dependance of e°, f° and g° 
(introduced in section [v} in the detector space (pi,p 2 ): 

e °M) = ^2 J e-^e^^^^ + qi)e- 9j ^dqidq 2 (58) 
Putting, gi = rcosO and 52 = rsinO we have 

e°(p,6» p ) = C* J o 27r J +o ° e - lr P cos ^- e ^rcos(r)e^drde 

C Jo + °° Jo* e- tr P cos ^- e ^de gi (r)dr (59) 



with C being the constant part. Now , taking 8 — 8 p — 8' the theta integral in eqn. (591 becomes 



27T p2ir-6 p 

e -irpcos(8-0 p ) _ / e -irpcos(9') ^gl 







e - lr pcos(9') del (6Q ) 

The last equation follows from the fact that the integral of a periodic function over its period is independent of the 
limits of integration. Thus e° is only a function of p: 

e°= ei (p) (61) 



Again, similar transformation for / yields, 



/•2?T Z' + OO 2 

f°(p,8 p ) = C / e- irpcos ^- e ^cos(8)sin(r)e^rdrd8 (62) 



jo 



Now, 



with C being the constant part. Using , 8 — 8 p = 8' the theta integral in (62 1 becomes, 

e - lrpC os(e') C0S ^ + e ^ de , = j^T-e P e -ir P cos(e') cos ( d ty os ( dp } de > 

~ f-T p 6p e- ir P co < e ">sin{d')sin{8 p )dd' 

2tt-0 p p2tt 

e- irpcos(e ' ) sin{8')d8' = e' lrpcos{0 ' ] sin{8')d8' . (63) 

-e p Jo 

Taking, cos{8') — x, the above eqn. becomes L 1 e~ lrpx dx + J", e~" px dx = 0. Therefore, 

f(p,e p ) = f x (p)cos{6 p ). (64) 
Proceeding exactly similarly one can show that 

9°(p,0 P )=fi(p)sin(6 v ) (65) 
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FIG. 4: The above figure represents the detector momentum plane. The positive pi and p 2 axes are represented by OD and 
OG respectively. OD is also the p[ axis. The p' 2 axis OC makes an angle 8 with OD. By DOE , we mean here the region of the 
plane bounded by OD and OE . 



2. Joint measurement probabilities 
Choosing the p[ , p' 2 axes according to fig. 4, the joint probability is given by , 



p(p'i>0,p' 2 >0) = p(pi,P2)dpidp 2 + p{puP2)dpidp 2 

Jpi,P2=0 J DOE 



(66) 



with, p(pi,P2) given by (|28j) . 

where the last integral represents integral on the region DOE in the figure 4. In polar coordinates the region DOE 
is given by, the set of points {(p, 9') : < p < oo, — — 9) < 9' < 0}. So, using eqn. (31 ) in (66) we have, 



p(p' 1 >0- / p' 2 >0)^(x\ 



—1+a — 
4 4 



Ix) 



p{p\,Pi)dp\dp 2 



DOE 



Now using the 9 dependance of e°, f° and g° from (611 , (64) and (65) in (28) we have, 
p(pi,P2)dpidp 2 = (x\(J™Jo { *~ 9) 



DOE 



ei(p?\ + \h(p?\cosH9) + |A(p) 2 | S m 2 (0))l s 
+2e 1 (p)f 1 (p)cos{9)a 1 + 2e x {p) h{p)sin(0)a 2 pdp<M)\x) 



Now, from eqn. (32) we have, 



/ (ej + ff)2irpdp = 1 
Jo 



We also have from the definition of a' (eqn. (34)) 



eihpdp = 



Hence 



p(p[>0,p 2 >Q) = ( X \ 



. 1 9 . 
( 2-^ 



— (1 + cos(9))a x + —sin(9)<T y 



(67) 



(68) 



(69) 



(70) 



(71) 
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and similarly, 



p(p[ < 0,p£ > Q) - ( X | 
p(pi>0,^<0) - (xl 
p(p'i<0,p 2 <0) - (xl 



2vr 

4" 



' "I 



(cos(8) 
(-cos(8) + l)cr 



Ix) 



-sin(8)(7 l 



Ix) 



2tt' 



^(1 +cos(6»))cr x - %m(6»)cr a \ : 



(72) 
(73) 
(74) 



At 8 = 5, we get back the POVM elements for the orthogonal case as derived in section |V| . For 8 i.e, joint 
measurement along two almost same directions, the probabilities p{p'i < 0,j?2 > 0) and >0,P2 ^ 0) vanish while 
p(p'i > 0,p 2 > 0) and p(p[ < 0,j»2 < 0) give the probabilities for single unsharp measurement along i, as expected. 

Remark: The joint measurement uncertainty relation is given by ([9| 



||3+6|| + ||a-6|| <2 

If is the angle between a and b , in the case a — b (like above) one has, 

1 

a < 



|swi(f)+cos(f)| 



(75) 



(76) 



The denominator in eqn. (76) is positive in 6 = [0,7r] with maximum y2 at 8 — tt/2. This is also our bound in the 
approximate joint measurement we have implemented, because for the symmetric inital detector state state case (see 
section VI A) approximate joint measurement happens in orthogonal directions as well. In fact we just showed that 
measuring momenta in any two directions yields approxiamte joint measurement of spin in those two directions. We 
saw earlier that in our case(see section VI) a is able to reach about 0.628 (see fig. la) . Thus, it is possible that one can 
have a different scheme with same but better quality of unsharp measurement in both the directions. For example, 
for 8 — it/4 the bound by the joint measurement inequality (76) is about 0.765. 



IX. SPIN DIRECTION FIDELITIES 

In this section we consider the spin direction fidelities, defined in ([12]) . We consider type 1 measurement as 
defined by the author in ref . |12) in which the pointer observables are taken to be the commuting components of some 
unit vector h. The measurement scheme considered by us in section [V] yields the momentum values p\ ,p2 (which 
can also be considered as components of p). Alternatively, in polar coordinates we can look at the magnitude of 
the momentum p and its angle 8 with the x-axis in the momentum plane. This angle uniquely fixes the direction of 
normalised momentum meter p, and this is what we take for h in this section. 

The joint unsharp measurement on the system through measurement on the commuting meter observables after 
the meters have interacted with the system, yields the POVM described by eqns. (33) and (36) .The effect of the 



measurement on the system ( a completely-positive map acting on the system density matrix) can be described by 
Krauss operators for the measurement T(p) . They are defined as follows by considering the average density matrix 
of the system pf after the measurement on the meters, considering all possible outcomes: 



Ps = / pdpd8T{p,8)\ X ){x\THp,8) 



(77) 



The probability of registering the outcome at an angle 8 in the pi , p 2 plane is given by p(8) = (x|£'(^)|x): with 
the angle POVM 



E{8)= / T^p,8)T(p,8)pdp. 
Jo 



(78) 



16 



The fidelities were verified to be expressable in terms of E and T as, (see 



Vi = in f\ x )eH sy M\ 

It = in f\x)eH sys (x\ 

Vd = inf\ x )eH, v M\ 
with T in our case given by, 



2tt 



(E{9)S.p + S.p(E(9))d9 



\x) 



T'n ■ STpdpdO 



\x) 



I ^(xl^SiTSi + SiT^SiT^dpde 



\x) 



T(p) = (p\e\^)l s + (p\m<j x + (p\g\^y 
= e°l s + f°a x +g a a y 



From (78) we have, 



r+oo 

E(0) = / ( 
Jo 



„0|2 , |i-0|2 



Hence, 



r+oo r+oo 

\f°\ 2 + \9°\ 2 )pdpls+ 2Re(e°f*)pdpa x + 2Re(e°g *)pdpa y 

Jo Jo 
a(6)l 8 + b(6)a x + c(9)a y ( say ). 



^(E(n)n.S + n.SE{h)) = ^[(b(0)cos{6) + c(6)sin(9))l s + a{9)cos(6)<T x + a{9)sin(9)a y ] 



So from eqn. (79 1, 

th = inf {x)eHsys d9[\{b{9)cos{9) + c{6)sin(9)) + \a{9)cos{9)( X \<j x \x) + \a{9)sin{0){ X \<J y \x)] 
= fa°° C e°f°*cos(9) + e° g°* szn(9)pdpd9 + J 27r (|e°| 2 + |/°| 2 + \g^)cos{e)pdpd9{ X \cr x \x) 

+ Io°° C(\e°\ 2 + l/T + \9°?)sm{9)pd P d6{ x \a v \x) 

Changing back to Cartesian coordinates we have, 

P2 



Vi 



e°f° \ 2 P1 2 ^ d Pl dp 2 

iPl+P2) 2 



e°g° 



00 (pi + pi) 2 



T dpidp 2 . 



(79) 
(80) 

(81) 



(82) 



(83) 



(84) 



(85) 



(86) 



The last two terms on the RHS of eqn. ( 85 ) vanish because the 3rd term is is odd in pi and the 4th one is odd in p 2 
Now using (70 1 we have, 



r)i = J e°fcos(9)pdpd9 + J e° g° sin(9)pdpd9 



/ eifipdp I cos 2 (9)d9+ I e 1 f 1 pdp 
J Je=o J 

2?r / eifipdp 



2tt 



=0 



\9)d9 



(87) 



Thus, 



m = (°°J 

Bound on a! : It was shown in [13] that ifo is bounded by the value s=.5 . From here we get the bound on a' to 
be about .64 which is almost exactly what is got in the simulation of a' (0.6292 in fig. la). The joint measurement 
uncertainty relation (43 1 allows a' to go till .707. This is because eqn. (43) refers to the most general approximate 
joint measurement in orthogonal directions without reference to any Arthur-Kelly kind of an implementation. Also 
eqn. (88) shows that the error of retrodiction (given in eqn. (21)) falls as the measurement becomes sharp. 
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A similar calculation shows that, 



Vf = inf {x)eHsys fpd9dp(x\[[co S (6)(\e°\ 2 + |/°| 2 - |. 9 °| 2 ) + 2Re(f°g°*) S in(6)]a x + [«n(0)(|e° 
+ |.g°| 2 ) + 2Re{f°g°*)cos(0)]<Ty + [e°f Q cos(9) + e°g sin{9)}l s }\x) 



\f° 



(89) 



and, 



Vd = in f\ x )eH sys / Pd9dp(x\ 



\f 



f\ 2 )h 



fO* 



o o* 
eg 



\x). 



(90) 



Thus, as in the case of original Arthur-Kelly model, the fidelities are independant of the initial sysytem statepQ. 
Again, because of the parity of various terms present there r] f tu rns out to be the same as rji 



From the probability normalisation condition given in eqn. ( 32 1 we have rjd = § , independant of initial apparatus 
state and rji ,r/f . 

This shows that any error-disturbance relationship between error of retrodiction/ prediction and error of disturbance 
does not hold for all choices of the direction pointer observable. 



X. APPROXIMATE JOINT MEASUREMENT FOR THREE QUBIT OBSERVABLES 

An approximate joint measurement observable for three unsharp qubit observables is defined by extending conditions 
given in ([!]) in the obvious way. Unlike the two observable case no necessary as well as sufficient condition for the 
approximate joint measurement on the marginal effect parameters in 72. 4 (like the ones derived for example in jlOj ) 
is known. No measurement uncertainty relation like equation ( 10 1 , but which is stronger than the relations applied 
pairwise is known. Here, we derive a necessary condition on the approximate joint measurement of three qubit 
observables which yields the familiar necessary-suffficient condition known in literature for the case of mutually 
orthogonal unbiased observables (to be discussed below). We show that it holds even when the observables are 
unbiased. Also when one of the measurements is trivial i.e the corresponding marginal effect is | representing 
equiprobable guesses of the values of the observable , our condition reproduces equation Q . 

The eight joint effects G +++ , G +4 , ,G have to satisfy six marginality conditions like 



Ti = G 



+++ 



G+A — + Gj L + Gj 



(91) 



A joint measurement of three observables gives rise to three two-observable joint measurements. Let G+ + denote the 
joint measurement marginal effect corresponding to the outcome (H — h) in directions 1 and 2. Then seven of the the 
effects (G +++ ,G + ^ , etc. ) can be written in terms of the three single-observable marginal effects , three pairwise 



joint measurement marginal effects and G+ 

G++- 
G-+- 

G__. 
G-+^ 

G___ 
G, 



as follows : 

G12 / - 



G 



= G 
= G 



+H 
23 
+- 

G w + 

- 

°++ 
23 
++ 
12 



G + + +; 

G+++ - 



G 4 



^<12 

- Lr + _ 

+ Gf_ 
-Gf_ 

- G+++, 

+ G u _ - G™ 

- Ti + G" 



, Y , 2 

1 +' 



G. 



++! 



G_i 



(92) 
(93) 
(94) 
(95) 
(96) 
(97) 
(98) 



It was shown in |10) that any set of effects for joint measurement of two observables (defined in eqn.|T])) can be written 
in the form 



G 



sgn(a)sgn(b) 



1 + ax + by + abZ + (abz + am + bn).a 



(99) 



l±Q+r 



and T± (y, n) 



0!) 



is true because for the two 



with a, b £ {1,-1} and the two single-qubit marginal effects given by , T±(x,fn) 
i±(y+ n .a) p sitfvity constraint being |a;| + m < 1 and \y\ + n < 1 . Condition 

observable case all the effects can be expressed in terms of one effect (say G ++ ) and the marginals (0) • The freedom 
in (Z, z) suffices to specify an arbitrary G++. 
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It then follows from equation [99] that the set of joint effects for the three observable case is of the form, 

G sg „(a)s a n{b)s a n(c){[Z], [z\) = {I + ax + by + cz + abZ x + bc:Z 2 + caZ 3 + abcZ 4 + {abz[ + bcz 2 + caz z 

+abczl + al + bfh + cn))/8 (100) 



with, a,b,c£ {1,-1}, the two-marginals G\ 2 = G(z\,Z\) , G23 
the one marginals being T±(x,l) = 1± ( x + Lcr '> ; T±(y, fh) = 1± (y+ m - rT 



G{z' 2 ,Z 2 ) , G13 = G(zs,Z 3 ) given by eqn 99 and 
T±(z, n) = 1± ( z +"- cr ) _ Given the positivity of 



the one marginals , the positivity of the eight joint effects places restrictions on the marginal effect parameters which 
are interpreted as measurement uncertainty relatons in contrast to the usual state uncertainty relations. 



A. Necessary condition 



From eqn. ( 100 1 , G +++ > gives, 

\\l + m + n + z*i+ £2 + Z3 + zl\\ < l + x + y + z + Z 1 + Z 2 + Z 3 + Z 4 



while G > gives, 



I - m — ft + z[ + z 2 + £3 — Z4 1 1 < 1 — x — y — z + Z\ + Z 2 + Z 3 



(101) 



(102) 



Equations (101) and (102) together can be interpreted as the collection of all points (z{ + z 2 + Z3) 6 1Z 3 , each 
of which lies within two spheres in 1Z 3 with their centres at (I + m + n + z\) and — (I + m + n + Z4) and radii 
(1 + x + y + z + Z\ + Z 2 + Z3 + Z4) , (1 — x — y — z + Z\ + Z 2 + Z3 — Z4) respectively. Thus, the distance betwee 
their centres should be less than or equal to the sum of their radii implying, 



I — in — n — zA\ < 1 + Z\ + Z 2 



(103) 



Similarly, by consideration of the other complementary pairs G pqr and C(_ p )(_ 9 )(_, r ) for p, q,r — ± we get three 
other equations, 



1 1 1 + rh — ft — z\ 1 1 < 1 + Z\ — Z% — Z% , 
1 1 1 — fh + n — z\ 1 1 < 1 — Z\ — Z% + Z3 , 
— I + rh + ft — z\ 1 1 < 1 — Z\ + Z% — Z3 . 



(104) 
(105) 
(106) 



B. Geometric intepretation 



The above inequalities ( 103 1 - ( 106 1 imply that we have a cuboid with edges 21 , 2m and 2ri with origin at the centre 
of this cuboid and solid spheres about points A (— I — rfi — n) , B(Z + m — ft), C (— I + rh + n) and D (I — fh + ?i)have 
to intersect so that the sum of the radii of the spheres is 4 (see fig. 5). That is we necessarily have: 



|| — I — fh — ft — ziW + \\l + fh — ft — 24II + ||/ — fh + n — za\\ + \\ — I + rh + n — < 4 (107) 

Now, for any given set of points in TZ 3 , the Fermat-Toricelli (FT) point of the set is defined to be the point which 
minimises the sum of distances from that point in 1Z 3 to the points of the set. The F-T point has been studied for 
quite long and its properties for any set of four non-coplanar points are well known. ([18J. For example, it is known 
that the F-T point either belongs to the set itself or else it is the point at which the gradient of the sum of distances 
vanish. Thus choosing 24 to be the F-T point of the set of points A,B,C,D one necesarily has (according to the 
definition of the F-T point), 

|| — I — fh — ft — z&f\\ + ||/ + fh — ft — z±p\\ + \\l — fh + ft — zipW + || — I + m + ft — zl F \\ < 4, (108) 
where z^p denotes the F-T point. 
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FIG. 5: Geometric interpretation of inequalities ( 104 1 - ( 106 1 . Solid spheres about the centres A (with OA = — (Z + fh + ft)) 
, B (with OB = (1 + fh - ft)), C (with 0~C = (-Z + fh + ft)) and D (with 0~D = (f- m + ft)) and radii (1 + Z\ + Z 2 + Z 3 ) , 
(1 + Z\ — Zi — Zz) , (1 — ^1 — ^2 + ^3), (1 — ^1 + ^2 — ^3) respectively have to intersect so that the sum of their radii is 4. 



C. Sufficiency 

In order to evaluate the F-T point for a particular set of points A,B,C,D, defined above we will use the following 
theorem from ref. |18j (Theorem 2.1), due to Lorentz Lindelf and Sturm. 

Theorem 2: Suppose zi is the F-T point for a set of points S n = {dl £ 1Z 3 : i = 1, 2, ...n}. Then either £4 belongs 
to the set S n or zX ^ S n . 

i) If zl £ S n then for z~i = Oj , for some j G {1, 2, ...n}, with || Y%=i(ijLj) pf^fjf II ^ 1 - 

ii) If Zi ^ S n then z~i is the point at which 1 1 iji^gT] 1 1 = 

Thus, for the first case the resultant of unit vectors to the FT point from other points of the set has magnitude less 
than 1. In the second case the unit vectors from the F-T point to the points of the set add up to the null vector. The 
condition for the second case is also the condition for the gradient of the function representing the sum of distances 
from z~i to dl, for i — 1,2, ...n, to vanish at Z4. 



1. I, in, n mutually orthogonal 



Suppose we are considering the case of approximate joint measurement in three orthogonal directions I , rfi and n. 
In this case we have, || — Z — m — n\ \ = \ \l + m — n\\ = \\l — m + ft\\ = || — I + m + n\\. Hence, at the origin(z4 = 0) of 
the cuboid (see fig. 5) we have, 



E 



(d* - z2) 



Zi\ 



= (— Z — m — n + Z — in + n + I + m — ri — Z + m + n)/(l 2 + m 2 + n 2 ) = 



(109) 



Hence equation ( 108 1 gives, 



< 1 



(110) 



|17j . The necessity of this condition was shown in (16] assuming 
= in eqn. (100)) by considering measurements by two parties on a 



This condition was shown to be suffficient in 
unbiasedness of the marginals (i.e x = y = z - 
singlet state. We have shown above that this is true for biased cases as well. Sufficiency of the above condition (jl 101) 



shows that condition ( 108 1 is sufficient for the case of approxiamte joint measurement in three orthogonal directions. 
Also note that condition ( 110 ) is stronger than pairwise conditions for two-observable joint measurement in orthogonal 



directions like (43) which when added produces the bound of 1.5 on the lhs of equation (110). 
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2. Reduction to two-observable inequality 

Suppose our joint measurement scheme is such that some approximate joint measurement on two observables is 
performed while the value of the third observable is guessed with probabilty of + being =4^ and that of - being ijp . 

This will correspond to n = i.e the corresponding marginal T±(z,n) = ^ 1= ^ z ^ . In this case the points A(— / — m) 
, B(7 + fh ), C(l — rh) and D (—1 + m) form a parallalogram of lenth |2Z| and \2m\ about the origin O . OA and OB 
lie opposite to each other and so does OC and OD. Thus, condition (ii) of the theorem 2 is satisfied and the origin is 
the F-T point. Hence, condition (108) reproduces equation ^ which shows that the bounds on the unsharpnesses is 



not more stringent than the two -observable case, as is to be expected. 

D. Arthur-Kelly model 

The case for three qubits proceeds exactly as that of two qubits with the Hamiltonian of the measurement interaction 



of the form, H = qi ® a x + q 2 ® a y + q$ eg) a z . The nature of the interaction is impulsive as in eqn. (24 1 and the 
corresponding unitary evolution is given by, 

U = e{qi,q 2 ,qz) <E) <r x + /(<?i, g 2 , <?3) ® v y + g(qi, ?2, fe) <8> <J Z (111) 

with, 



e(9i,92,?3) = cos{^J ql + ql + ql) (112) 

f[qi,Q2,Q3) = m — n \ 2 =-, (113) 

, \ ■ sin{yj q\ + ql + ql) 

5(91,92,93) = — 2 2 2 =^, (H4) 

1 / s sin{ytql + ql+ql) 

h(qi,q2,q3) = m 7= „ 9 =-■ (115) 

+ «| + 9a) 

The intial joint state of the three detectors and system is \ipi n ) = ® ^2) ® |"03) ® lx)- As before, the initial 
detector states are gaussians: if>i(qi) = i^ e '^)^ , ^2(92) = i^ e ~ 2%)^ and ^3(93) = (^ e ~l^)'- 

The detector outcome ((pi > 0,p 2 > : pz > ) is taken here to correspond to the outcome (+,+,+) for the joint 
unsharp measurement of the system observables a x , <r y , a z ; (pi > 0,p 2 > 0,ps < 0) to (+,+,-) and so on. The 
POVM elemnts corresponding to the outcomes (+++), (++-) , etc. of the joint unsharp measurement of a x , a y 
and a z have a similar structure to that for the case of two observables : (+ + +) O g(7 + + \ a y + % a z) , 

(+ + -) U T + + b s a v - » ( ) ° W - - b 8°v - with ' 

a'= / / 4(f e )d Pl dp 2 dp 3 , (116) 

J pi — J P2——00 J p 3 — — 00 

&' = / / / 4( 5 °e°)dp 1 dp 2 dp3, (117) 

J p 2 — ^ Pi — — OO J P3 — — OO 
/»+oo /*+00 /"+OO 

c' = / / / i(h°e°)dpidp 2 dp 3 . (118) 

J PS—O J P1 — — OO J P2 — — OO 



-Z. Symmetric case with a=b=c 

With the initial joint pointer state a symmetric Gaussian, as in the two observable case, we have for the unsharp 
measurement along any three directions h,rh,l, the marginal effects T™ = ~(J + a'n.a) , T" 1 = ~(J + a'rh.a) and 
T' + = i(J + a'Z.c?) . Thus from eqn. 

scheme, a' is seen to be able to reach upto about 0.49 . It cannot reach the bound 0.577, as for two-observable joint 
measurement. (see fig. 6). 



110) we have a' to be bounded by — 0.577 . Numerically , for our 
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I n« - 



FIG. 6: A plot of a' vrs a for symmetric Gaussian initial state 



For approximate joint measurement in the directions I = (1,0,0) , rh = (cos(0), sin((f>), 0) , n = 
(sin(6)cos((j)i), sin(9)sin((j)i), cos(9)) with 9 = 0.414-7T, <\>\ = 0.1597T and cj) = 0.4777T, the Fermat- Toric elli point of 
the set of points A,B,C,D, as defined before, is seen to be at / + m — n . This yields, from inequality ( 108 ), a' < 0.667. 
Thus, as for the two observable case, more freedom is available for unbiased measurement in non-orthogonal directions. 



XI. CONCLUSIONS 



Non-commuting observables cannot be measured joinly. However it is possible in the POVM formalism to do 
joint measurements of (generally non-commuting) observables which are approximations of the actual non-commuting 
observables and called unsharp observables. We have considered in this paper approximate joint measurement of two 
and three qubit observables separtely, through an Arthur-Kelly like model for qubit observables. This model comes 
naturally when one considers a Stern-Gerlach setup with a linear magnetic field. In Stern-Gerlach setup the momenta 
of the atoms acts as pointer observables for their spin degrees of freedom. Considering approximate joint measurement 
of a x and a y through this model, we have shown here numerically that measurement uncertainty relations derived 
elsewhere(see ref.[7j) hold. It has also been shown how increasing the relative sharpness of the initial momentum 
wavefunctions of the detectors leads measurement of one observable to become almost sharp while that of the other 
to become almost trivial. Effect of initial detector states on the post-measurement state of the sysytem has also been 
considerd. 

We have considered two different characterisations of unsharpncss. Firstly, by comparing the probability distribution 
of the observable to be approximately measured with that of the approximate observable. Secondly, by considering 
the alignment of the momcmtum direction with the spin observable in the Heisenberg picture. We have shown that 
for the case in which both the observables are approximated equally well, the corresponding measures of unsharpness 
are proportional. For our choice of pointer observable the measures checking the alignment and disturbance due to 
measurement do not seem to satisfy an error-disturbance relationship. 

We have expounded the connecion between the symmetries of the underlying Hamiltonian for measurement inter- 
action and initial detector states with the joint measurement in a lemma. This was first stated by Martens et al.([3]). 
It has then been used to perform approximate joint measurement in arbitrary directions. The POVM elements cal- 
culated for the same match with that which were found earlier in the orthogonal case. They also turn expectedly to 
that of a single unsharp measurement when the directions are taken to be almost same. 

For the case of joint unsharp measurement of three qubit observables we have given a necessary condition to be 
satisfied by the parameters of the marginal POVM elements . This condition has been derived from certain geometric 
considerations involving the so called Fermat-Toricelli point . The condition is sufficient for the case of three orthogonal 
observables and identical to the only necessary-sufficient condition known for three -observable joint measurements. 
Our proof shows that this also holds for biased unsharp measurements , namely those in which the probabilty of 
obtaining 'up' or 'down' is different for a maximally mixed state. 

The measurement scheme employed by us for joint measurement in non-orthogonal directions cannot take advantage 
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of the greater freedom available for better approximation compared to the orthogonal case. It will be interesting to 
see such a scheme in the Arthur-Kelly setup that is able to get close to the bound set by (eqn{9| for arbitrary 
directions. We have shown that the retrodictive fidelity r\i and the unsharpness a! are proportional for the symmetric 
joint unsharp measurement case that we have considered. It may be possible to connect the two pictures in a 
more general settig starting with certain symmetries of the Hamiltonian and initial detector states. The problem of 
determining necessary-sufficient conditions for the most general joint measurement of three observables by extension 
of our approach or otherwise is also left open. 
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